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1.  Introduction 


There  is  a  longstanding  controversy  in  the  ballistic  community,  largely  undocumented, 
regarding  the  correct  value  of  the  internal  energy  of  formation  of  condensed  phases  as  it  is  used 
in  thermodynamic  equilibrium  calculations.  Recently  this  problem  again  surfaced  at  a  meeting  of 
NATO  AC/225  (Panel  IV/SD.2)WP9  which  was  working  on  the  standards  agreement,  STANAG 
4400,  "Derivation  of  Thermochemical  Values  for  Interior  Ballistic  Calculations"  (France,  September 
1990).  At  this  meeting  the  results  of  thermodynamic  equilibrium  calculations  using  two  differing 
interpretations  (Freedman  1988;  Volk  and  Bathelt  1990)  and  corresponding  values  of  the  internal 
energy  of  formation  were  presented.  This  particular  problem  does  not  arise  for  constant  pressure 
calculations  where  the  relevant  thermodynamic  quantity,  the  enthalpy  of  formation,  is  tabulated  and 
can  be  used  directly.  Equilibrium  calculations  in  which  the  reactant  is  a  gun  propellant,  however, 
are  usually  performed  for  a  constant  volume  system,  since  this  situation  is  assumed  to  better 
simulate  the  conditions  in  a  gun.  The  internal  energy  of  the  propellant  is  then  the  quantity  that  is 
usually  required;  it  is  presumed  necessary,  therefore,  to  calculate  the  internal  energy  of  formation 
from  the  tabulated  value  of  the  enthalpy  of  formation,  and  it  is  here  that  the  problem  is 
encountered.  The  purpose  of  this  report  is  to  clarify  the  meaning  of  the  internal  energy  of 
formation  as  it  relates  to  thermodynamic  equilibrium  calculations,  especially  as  it  pertains  to 
calculations  which  employ  the  standard  computer  programs  now  in  use. 


2.  Background 

Two  examples  from  the  literature  are  selected  below  to  illustrate  the  divergent  views. 
Before  proceeding  with  these  examples,  some  basic  definitions  and  conventions  are  first  reviewed. 

The  enthalpy,  H,  is  a  convenience  function  applicable  to  constant  pressure  processes  and 
is  defined  as 

HsU  +  PV  (1) 

where  U  is  the  internal  energy,  P  is  the  pressure,  and  V  is  the  volume.  From  this  it  follows  that 

AH  =AU  +  A(PV) 
and 

AHy  =  AUy  +  A(PV) 

where  AHy  and  AUy  are  the  enthalpy  and  internal  energy  of  formation  respectively. 
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The  energy^  that  is  customarily  tabulated  is  the  standard  enthalpy  of  formation  (sometimes  called 
the  heat  of  formation),  AHy°,  defined  (Chase  et  al.  1985)  as  "...  the  increment  in  enthalpy 
associated  with  the  reaction  of  forming  the  given  compound  from  its  elements  in  their  reference 
states,  with  each  substance  in  its  thermodynamic  standard  state  at  the  given  temperature."  By  this 
convention,  for  all  elements  in  their  standard  state,  AHy°=0  at  all  temperatures.  The  superscript 
denotes  that  the  pressure  is  for  the  standard  state.  In  older  tabulations,  the  standard  state  (pressure) 
is  usually  1  atm;  more  recent  compilations,  e.g.  the  JANAF  tables  (Chase  et  al.  1985),  use  a 
pressure  of  1  bar  (=  0. 1  MPa)  for  the  standard  state,  and  convert  the  older  values  based  on  1  atm 
to  this  new  standprd  state  pressure.  It  is  also  customary  and  convenient  to  designate  an  "absolute" 
energy  scale  by  assigning  the  enthalpy  of  the  stable  form  of  elements  at  the  standard  pressure, 
hereafter  designated  as  p°,  and  a  reference  temperature,  usually  298.15  K,  a  value  of  zero,  /.e., 
H°298.i5  =  H(298.15  K,  p°)  =  0.  There  actually  is  no  absolute  energy  scale;  however,  by 
designating  the  standard  state  (p°)  and  the  reference  temperature,  and  assigning  the  enthalpies  of 
the  reference  elements,  an  effective  "absolute"  scale  is  created.  Hereafter  in  this  report,  the 
designation  assigned^^  applied  to  either  an  enthalpy  or  internal  energy  will  mean  an  energy  which 
is  based  on  this  chosen  reference  state,  i.e.,  H(298. 15  K,  p°)=0  for  elements  in  their  stable  form 
of  aggregation  at  298.15  K  and  p°.  For  a  change  of  state,  a  negative  value  of  AU  or  AH  means 
that  energy  is  evolved  by  the  process. 

The  relationship  between  the  standard  enthalpy  of  formation,  as  it  is  defined,  and  the 
assigned  enthalpy  is  illustrated  in  Figure  1.  (AHy°)j  is  the  enthalpy  of  formation  at  some 
temperature  T  for  which  the  reactant  elements  of  the  formation  reaction  have  some  arbitrary 
enthalpy  value.  Also  shown  is  the  same  (AH^°)j  on  the  assigned  scale,  where  it  is  displaced  by 
an  amount  equal  to  (H°^^^,^„,^>r-(H°„^,^„„)298.i5K  =  <-^reactants>T-  The  quantity  (AH/)t  has  the 
same  value  on  any  energy  scale;  in  general,  however,  the  assigned  enthalpy  H°j  ^  (AHj°)j  for 
a  temperature  other  than  the  reference  temperature,  in  this  case  T  =  298.15  K,  of  the  assigned  scale. 
[The  situation  illustrated  in  Figure  1  is  typical  for  (AHy°)j  when  T>298.15K].  Also  note  that, 
owing  to  its  definition,  discontinuities  can  occur  in  AHy°  if  the  stable  form  of  aggregation  of  an 
element  passes  a  phase-transition  temperature  (Chase  et  al.  1985) . 

U  and  H  can  be  expressed  as  functions  of  any  two  other  state  properties.  Choosing  the 
appropriate  natural  variables  gives 


and 


U=U(T,V) 

(2a) 

H=H(T,P). 

(2b) 

For  ideal  gases,  H  and  U  are  rigorously  only  functions  of  T.  In  general,  however,  as  indicated  in 
Eq.  2,  H  and  U  also  have  a  pressure  and  volume  dependence,  respectively;  if  the  initial  conditions 


"  The  term  energy  (italicized)  will  be  used  in  this  report  to  refer  to  either  internal  energy  or  enthalpy. 
This  designation  is  also  used  in  the  NASA-Lewis  documentation  (Gordon  and  McBride  1971). 
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Assigned  Scale 


Figure  1.  Schematic  energy  diagram  showing  the  relationship  of  the  enthalpy  of  formation  at  some 
arbitrary  temperature.  T.  to  the  assigned  energy  scale.  The  enthalpy  of  formation  at  298. 15K  is 
also  shown  for  comparison. 


for  a  calculation  are  other  than  those  of  the  tabulated  standard  enthalpy  of  formation,  then  its  yalue 
or  the  corresponding  internal  energy  must  be  adjusted,  as  necessary,  for  the  speciiled  initial 
conditions  of  temperature,  pressure,  or  yolume.  For  example,  on  the  assigned  energy  scale,  for 
standard  pressure  p°,  and  reference  temperature 


where  Cp°(T)  is  the  standard  state  heat  capacity.  Howeyer,  for  P  p°, 


dP 

T 
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Unless  otherwise  noted,  in  the  following  discussions  all  compounds  will  be  considered  to 
be  in  the  standard  state,  i.e.,  P=p°,  and  at  the  reference  temperature  T=298.15  K.  The  two 
representative  approaches  are  now  summarized.  Both  cases  have  the  same  expressed  purpose: 
deduce  the  value  of  the  internal  energy  of  formation  from  the  enthalpy  of  formation  for  a 
condensed  phase. 

2.1  Case  1 

In  this  approach  (Freedman  1988),  it  is  noted  that  the  internal  energy,  U,  is  the  appropriate 
variable  for  constant  volume  calculations,  and  that  the  internal  energy  of  formation,  AUy,  is  the 
"correct  input  datum",  it  is  pointed  out  that  the  difference  between  the  and  AUy  i.e.  the  PV 
contribution,  is  small,  and  is  less  than  120  cal/mol  for  typical  propellants;  this  value  is  calculated 
for  a  pressure  of  1  atm  and  a  maximum  value  of  0.005  m^  for  the  molar  volume  of  a  propellant 
ingredient.^  Furthermore,  it  is  argued,  this  small  a  difference  is  less  than  the  uncertainties  (ca.  0.3- 
2  kcal/mol)  of  tabulated  values  of  the  enthalpies  of  formation.  For  condensed  phases,  therefore, 
it  is  correct  to  use  the  value  of  the  enthalpy  of  formation  for  the  internal  energy  of  formation  when 
the  latter  quantity  is  required. 

2.2  Case  2  ^ 

In  the  other  approach  (Volk  and  Bathelt  1990),  the  need  for  the  energies  of  formation,  rather 
than  the  enthalpies  of  formation,  when  performing  constant  volume  calculations  for  gun  propellants 
is  also  noted.  The  energy  of  formation,  AU^,  is  calculated  using  the  standard  formula: 

AUy  =  AH^  -  An  R  T  =  AH^  -  An  (0.59248)  (kcal/mol), 

where  An  is  the  difference  in  the  number  of  moles  of  gases,  i.e.  ^products  *  ^reactants'  ^ 
universal  gas  constant,  and  T  is  the  absolute  temperature.  [RT=0.59248  kcal/mol  for  T=298.15K]. 
The  term  AnRT  is  a  standard  approximation  for  A(PV)  assuming,  as  in  the  first  case,  that  the  PV 
product  of  condensed  phases  is  negligible.  A  sample  calculation  is  performed  using  this  approach 
for  the  internal  energy  of  formation  of  dibutyladipate,  a  liquid,  at  T=298.15  K,  chemical  formula 
C14H26O4,  AHy=-181.7  kcal/mol,  for  the  formation  reaction 

14  C  (solid)  +  13  H2  (gas)  +  2  O2  (gas)  CJ4H26O4  (liquid). 


Since  An  =  0  -  15  =  -15,  the  AU^  is  readily  calculated  to  be  -172.8  kcal/mol.  The  difference,  in 
this  approach,  between  AU^  and  AH^  is  AnRT  =  8887.6  cal/mol. 


To  take  an  actual  example  for  comparison,  at  P=1  atm.  for  HMX  which  has  a  density  p=1.9  g/cm^  a  molecular  weight 
MW=296.17  g/mol,  and  using  V=MW/p,  the  PV  product  is  just  3.78  cal/mol. 
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The  two  methods  for  calculating  Ally®  above  give  results  for  the  related  quantity  enthalpy 
of  formation  -  internal  energy  of  formation  which,  when  compared,  are  seemingly  in  sharp 
disagreement  Case  1  estimates  a  maximum  difference  of  120  cal/mol  for  a  typical  propellant  The 
method  used  in  Case  2  gives,  for  a  liquid,  8887.6  cal/mol.  Since  both  methods  are  dealing  with 
a  condensed  phase  substance,  this  large  a  disagreement  120  cal/mol  versus  8887.6  cal/mol,  is  not 
expected. 

3.  Discussion 

The  resolution  of  this  discrepancy  can  best  be  addressed  by  reference  to  a  schematic  energy 
level  diagram.  For  convenience,  and  by  way  of  illustration,  the  internal  energy  and  enthalpy  of 
dibutyladipate,  which  was  calculated  in  the  second  example  above,  will  be  used.  In  Figure  2  the 


C(s)  ©2 


I14C(s)  +  13H2(fl)  +  202(flf)] 


U=-0.128/moh^ 
U=-593/mol  _ 


^2  ^2 


U=-8889.4 


(solids) 

(gases) 


H=0 


Reactant 


CO 

o 


U=H-PV{//q) 


Product 

Cl4«2604W 


H=-181700/mol  | 


Figure  2.  Schematic  energy  diagram  showing  the  enthalpy  and  internal  energy  states  involved 
in  the  formation  dibutyladipate. 


energies  of  the  various  compounds  are  drawn  schematically  on  a  diagram  that  is  only  intended  to 
convey  large  and  small  differences  in  the  quantities,  rather  than  provide  a  precise  scale.  The 
elemental  species  needed  to  form  the  dibutyladipate  are  listed  at  the  top  of  the  figure  at  the  zero 
energy  reference  state.  In  addition,  the  reactant  composition  for  the  formation  reaction  is  also 
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indicated,  and  is  included  in  brackets  for  identification;  it  also  has  an  H=0,  which  is  just  the  sum 
of  its  component  enthalpies.  Just  below  this  level,  at  -0.128  cal/mol,  is  a  representative  U  for  an 
elemental  solid,  where  the  particular  value  chosen  is  for  carbon.  This  and  subsequent  U  values  are 
all  calculated  using  Eq.  1.  The  small  difference  between  the  internal  energy  and  enthalpy  of  the 
solid  carbon  is  just  the  PV  product  term;  this  was  calculated  using  (Lange’s  Handbook  of 
Chemistry  1979)  a  density,  p=2.267  g/cm^,  a  molecular  weight  MW=12.01 1  g/mol,  and  the  relation 
V=MW/p.  The  third  level  down,  U=-593  cal/mol,  corresponds  to  the  PV(=RT)  of  one  mole  of  an 
ideal  gas.  The  next  (fourth)  level  at  U=-8889.4  cal  is  the  total  internal  energy  of  the  reactants 
needed  to  form  dibutyladipate.  Since,  as  noted  above,  H=0  for  this  reactant  composition,  U  is  just 
the  negative  of  the  sum  of  the  PV  product  of  each  species,  where  each  mole  of  (assumed  ideal)  gas 
contributes  RT=593  calories  and  each  mole  of  solid  carbon  contributes  0.128  calorie.  These  four 
levels  schematically  represent  all  the  internal  energy  and  enthalpy  states  of  the  reactants. 

The  next  two  levels  represent  the  enthalpy  and  internal  energy  for  the  product  compound 
dibutyladipate.  The  upper  level  of  the  two,  at  H=-181.7  kcal/mol,  is  derived  from  the  tabulated 
value^  of  the  enthalpy  of  formation  and  the  change  of  state  given  in  Case  2  above,  where,  for  the 
dibutyladipate,  H  =  AHy  +  H^eactants  -  -181.7  +  0.0  =  -181.7  kcal/mol.  The  next  slightly  lower 
level  represents  the  corresponding  internal  energy  value  of  dibutyladipate.  For  this  level  U=H- 
PV(//^),  where  'W(liq)  is  the  molar  volume  of  the  liquid  phase  dibutyladipate.  As  implied  in  both 
cases  presented  above,  at  pressure  p°,  the  value  of  PV(liq)  is  expected  to  be  much  less  than 
^product  -  -181700  cal/mol.  Neglecting  the  small  PV(/i^)  contribution  gives  U=H. 

The  levels  in  Figure  2  represent  the  reactant  and  product  assigned  energy  states  for  the 
formation  of  dibutyladipate.  These  levels  are  duplicated  in  Figure  3.  Additionally  in  Figure  3, 
vertical  arrows  are  drawn  to  indicate  the  states  involved  in  forming  the  AUy  and  AH^  for 
dibutyladipate.  Referring  to  Figure  3,  and  using  the  states  involved  in  the  formation  reaction  of 
Case  2,  it  follows  that,  for  the  product,  H=AHy^  exactly,  however,  U=AUy-(PV)^^^,^,.  Because 
of  the  choice  of  reference  state,  the  assigned  value  of  the  product  enthalpy,  Uproducr 
enthalpy  of  formation,  AHy,  are  always  equal.  In  general,  however,  the  assigned  internal  energy, 
U product  internal  energy  of  formation,  AU_^,  are  not  equal,  differing  by  iP^)reactanp  which 

is  not  usually  negligible  if  gaseous  reactants  are  involved. 

In  summary,  for  compounds  which  are  used  as  a  reactant  component  in  an  equilibrium 
calculation,  at  the  standard  state  pressure  and  reference  temperature  (here  298. 15K),  on  the 
assigned  energy  scale,  H°  =  AHj°  and  U®  ^  AUy°;  also,  for  condensed  phases  only,  usually, 
U°  =  H°. 

Two  commonly  used  thermodynamic  chemical  equilibrium  codes,  NASA-Lewis^  and 
BLAKE  (Freedman  1981),  calculate  a  standard  enthalpy,  denoted  H®(T),  which  is  a  function  of  the 
temperature  and  a  set  of  coefficients.  These  H®(T)  values  are  assigned  enthalpies  from  which 
assigned  internal  energies  are  calculated.  These  codes,  therefore,  use  an  assigned  energy  scale  for 
referencing  both  the  enthalpy  and  internal  energy  of  reactants  and  products;  the  input  data,  i.e.  the 
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internal  energy  or  enthalpy  of  the  propellant  components,  must  correspond  to  this  scale.  Formally, 

U  and  H  represent  these  assigned  energies  while  AU  and  AH  do  not 

4.  Conclusions 

1.  Because  the  choice  of  reference  state  is,  by  convention,  based  on  the  reactants  of  the 
formation  reaction  having  an  assigned  enthalpy  value  of  zero,  it  follows  that  at  the 
reference  temperature,  for  the  product  of  formation,  H°  =  AHj°. 

2.  The  quantity  (AHy®)j  has  the  same  value  on  any  energy  scale;  in  general,  however,  the 
assigned  enthalpy  H°j  *  (AHy°>j-  for  a  temperature  other  than  the  reference  temperature 

3.  Because  of  this  assigned  energy  scale,  in  general,  for  the  same  product  also  at  the  reference 
temperature,  U°  *  AU^°. 

4.  The  approach  in  Case  2  above  correctly  calculates  AUy°,  but  this  is  not  the  quantity  that 
is  likely  to  be  required  by  a  ttermodyntunic  equilibrium  code;  for  instance,  it  is  not  the 
quantity  required  as  input  by  two  widely  used  codes,  BLAKE  and  NASA-Lewis. 
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5.  The  relationship  between  the  internal  energy  and  enthalpy  of  condensed  phases  as  stated 
in  Case  1  is  also  correct,  however,  it  pertains  to  the  assigned  energies,  U°  and  H°,  and  not 
to  AU/  and  AH/. 

6.  In  general,  for  condensed  phases  referenced  to  the  conventional  energy  scale  having  H°=0 
for  elements  in  their  standard  state,  U°  =  H°  and  AU/  *  AH/. 
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